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A new type of boundary dynamics is proposed to describe the interface that sweeps space to 
collect distributed material. Based upon geometrical consideration on a simple physical process 
representing a certain experiment, the dynamics is formulated as the small diffusion limit of MuUins- 
Sekerka problem of crystal growth. It is demonstrated that a steadily extending finger solution exists 
for a finite range of propagation speed, but numerical simulations suggest they are unstable and the 
interface shows a complex time development. 
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I. INTRODUCTION 

The boundary dynamics has been providing interest- 
ing problems for physics and mathematics. A famous 
one is the problem of crystal growth from a supercooled 
melt Pi ; As the melt is solidified and a crystal grows, a 
flat interface between the solid and the liquid phases be- 
comes unstable due to the coupling of the solidification 
process with the diffusion of latent heat generated at the 
interface (MuUins-Sekerka instability^). This results in 
the fascinating variety of dendritic growths of crystal un- 
der the interplay with anisotropic surface tension. An- 
other example is the viscous fingering, which appears 
when the air is injected into a viscous fluid (Saffman- 
Taylor instability '3', 1). The viscous fluid is displaced 
by the pressure gradient, and the pressure field is gov- 
erned by Laplace equation with proper boundary condi- 
tions. Since Laplace equation is the diffusion equation 
with the infinite diffusion constant, the viscous fingering 
is the large diffusion limit of the crystal growth. 

A new example of boundary dynamics is presented 
by a simple experiment of Yamazaki and MizuguchiQ; 
The mixture of water and corn starch powder is sand- 
wiched between two glass plates. After a several hours, a 
labyrinthine pattern of dried corn starch will be formed 
when the water is evaporated from the gap of the glass 
plates. The system is two-dimensional, and the pattern 
is formed by the water-air interface line as it sweeps the 
system to collect the granules along it by means of the 
surface tension. The grains are simply accumulated along 
the interface line to give friction against the interface mo- 
tion, and eventually get stuck with the glass plates. 

In the above experiment, the granules play an analo- 
gous role to that of the latent heat in the crystal growth; 
They are distributed in the wet (or melted) region, but 
show up when the interface passes; The interface speed 
is controlled by the granule density (or temperature) at 
the boundary. The difference is that the granules do not 
diffuse while the heat does. In this sense, this sweeping 
dynamics is the small diffusion limit of MuUins-Sekerka 
problem; That is the opposite limit to Saffman- Taylor 
problem but has not been investigated yet in detail. 

Analogous instabilities exist for these phenomena; A 
protruded part of interface advances faster because the 



generated heat diffuses faster (MuUins-Sekerka), the pres- 
sure gradient is larger ( Saffman- Taylor ) , or the accumu- 
lated granules are diluted over the elongated interface at 
the convex region (sweeping dynamics). 

Similarities is seen in the phase field model proposed 
for this sweeping phenomenon j^, Q- The model con- 
sists of two fields: the phase field and the coupling field; 
The coupling field represents the granular density, in- 
stead of the temperature in the crystal growth[3. 19I [lol| . 
It has been demonstrated the model is capable to re- 
produce some feature of the patterns obtained in the 
experiment 0- 

In this paper, we will construct the model of the bound- 
ary dynamics for the sweeping interface based on geomet- 
rical and physical considerations, and study its behavior. 



II. SYSTEM CONFIGURATION AND 
CO-ORDINATE 

Let us start by defining the Cartesian coordinate in the 
two dimensional system near the interface between the 
swept (dry) and the unswept (wet) region as in Fig^a). 
In the unswept region, the granules are distributed at the 
area density p, which we assume constant in this paper, 
for simplicity. 

Suppose that the interface position at the time t is 
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FIG. 1: Schematic diagram of the interface between the swept 
and the unswept regions, (a) Coordinate system with the 
normal and tangential vectors, (b) Steadily extending finger. 
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represented parametrically by 

r{t,s) = {x{t,s),y{t,s)) 



(1) 



with the parameter s. Let £ be the natural coordinate 
along the interface, then the length element d£ is 



de = \J{dxldsf + {dy/dsf ds . 



(2) 



We define the tangential and the normal unit vectors of 
the interface, t and n, by 



i = dr/dl , n = £ X t , 



(3) 



where z denotes the unit vector along the z axis perpen- 
dicular to the system; The normal vector n is pointing 
into the unswept region. The curvature k, is defined as 



dt „ 



(4) 



which is positive when the interface is convex toward the 
unswept side. 

In the drying process, the interface moves upwards. 
Using an appropriate parameterization of s in Eq.(^, 
dr/dt can be made parallel to n, then the normal speed 
of the interface motion w„ is defined by 



dr{t, s) 
dt 



u„ it,s)n{t, s). 



(5) 



III. DYNAMICS OF SWEEPING INTERFACE 

We now consider the sweeping dynamics where the 
granules are swept by the interface; They are accumu- 
lated along the interface and conveyed by it in its normal 
direction. If we ignore the width of the region where the 
granules are accumulated, then the accumulated granules 
are described by the line density a along the interface. 

In the case where the granules are simply accumulated 
along the interface and do not diffuse at all, the equation 
for a is determined geometrically and should be given by 



da{t, s) 
dt 



Vn{t,s) p - K{t, S)a{t, S) 



(6) 



The first term of the right hand side simply represents 
the sweeping accumulation. The second term comes from 
the change of the interface length as it advances; The 
length element along the interface increases by the factor 
(1 -|- Kv„A.t) during the short time period At, thus the 
line density decreases by the factor of its inverse. 

The interface speed u™, on the other hand, is given by 
the product of the following two factors, i.e. the driving 
force to the interface, and the mobility of the interface; 
(i) The driving force comes from the pressure difference 
AP between the wet and the dry regions; As the water 
evaporates, the volume of the wet region tends to shrink. 



then the interface recedes due to the pressure difference. 
When the interface is curved, the driving force is given by 
the effective pressure difference (AP — 7k) with 7 being 
the surface tension, thus is proportional to the factor 
(1 — qk) with the capillary length a = 7/AP. (ii) The 
mobility depends upon the granule density a and should 
be a decreasing function of it because the granule friction 
with the glass plates resists the interface motion. The 
mobility becomes zero at dst when the interface gets stuck 
with granules. 

In the simplest case where no other scales of a are 
involved, we can write down the equation for the interface 
speed as 



Vn = 1'o/(o-/(Tst)(l - an) 



(7) 



with the characteristic speed vq and the dimensionless 
mobility f{x), which is a decreasing function with /(O) — 
1 and /(I) = 0. 

Eqs.(|HJl and {Tj) define the boundary dynamics of the 
sweeping interface, which shows morphological instability 
analogous to the crystal growth dynamics; The part of 
the interface with k > tends to advance faster when 
/' < 0. 

The model contains four parameters, p, vq, (Jst, and a, 
from which we can define the stuck-in distance 1st and 
the stuck-in time tst by 



4t = O-St/P, tst = £st/vo, 



(8) 



which are the distance and the time that the flat interface 
can advances before it gets stuck. 

In order to simplify the expressions, in the following, 
we employ the dimensionless unit system where £st = 
tst = cTst = 1, then Eqs.® and {Tj) are in the form of 



da{t, s) 
dt 



u„(t, s) 1 - K{t,s)a{t,s) 



Vnit,s) = f(^a{t,s)^(l^ Rn{t,s) 
with the only one dimensionless parameter 
R = a/£st = ap/ast- 



(9) 
(10) 

(11) 



This is the ratio of the capillary length to the stuck-in 
distance and is a measure of the effect of the surface 
tension. 

The functional form of the dimensionless mobility /(cr) 
in Eq. 1)10(1 should reflect the physical mechanism how the 
interface slows down and gets stuck due to the accumu- 
lated granules. Here, we employ the simplest form 



(1 - cr) for < (T < 1 
for (T > 1 



(12) 



IV. SIMPLE SOLUTIONS 



Now, we study the interface dynamics based on 
Eqs.®,lTni), and 
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flat interface 



when 



For the flat interface, k — 0, thus the solution is easily 
obtained as 

Vn{t) = e-' = Vfit), cr(t) = l-e^* = af(t). (13) 

The flat interface can advance only by £st(=l) before it 
gets stuck because it simply accumulates material. 



B. steadily extending flnger 

A possible mode of steady advance is the extending 
one-dimensional finger of the width of 2£st by shoving 
the granules asides(Figmb)). 

This type of steadily extending finger solution can be 
obtained as follows. Suppose the finger is extending in 
the y direction with the speed V, then Eqs.lPJ and ifTIHl 
become 

-Vcr'y' = {1-k(t)Vx', (14) 
Vx' = (1 -cr)(l -i?K), (15) 

where the primes denote the derivative by the natural 
coordinate £. Then, x' and y' are related as 



1, 



(16) 



thus, Eqs. 1)14(1 - can be solved for x, y, and a as 
functions of £ under the physical boundary conditions 



{x,y) 



ix,y) 
(Tl,-oo), 



= (0,0) 

cr = 1 



at 
at 





±oo. 



for a given set of R and V. 

For a given R, a steady solution is possible for a finite 
range of V; The condition at the tip that x' = 1 and 
y' = at ^ = in Eqs.ldll and leads to the allowed 
range of V, 



V < {1-VRf. 



(17) 



Eas. (|14() - (|16l) can be solved numerically. Fig[21 shows 
some of the steady solutions for the finger shape and the 
line density cr for i? = 0.2 with V = 0.1, 0.2, and 0.3. 
One can see that the finger tip is sharper for the larger 
extending speed V. 

The line density a shows rather intriguing behavior as 
a function of £; The tip density is smaller for the larger 
V, but for V larger than a certain value for a given R, 
a becomes singular at the tip and eventually develops a 
cusp. Actually, by expanding a around the tip, it can be 
shown that tr becomes non-analytic 



with 



If, = 



a[£) ~ ct(0) cx £'f 



2^(1- V" + i?)2-4i? 



[\-V -R)- + 



4i? 



(18) 



(19) 



1 

^+3 



hR - 4V-R2 + 3i?j < y < (1 - \fRf 



(20) 



The singularity at the tip in the line density is a peculiar 
result of the model without diffusion. 



V. ULTRA-VIOLET CATASTROPHE IN THE 
DIFFUSIONLESS MODEL 

We have derived the steady finger solutions, which are 
smooth and analytic except at the finger tip, but in a 
general time development, the diffusionless feature of the 
model seems to cause the ultra-violet catastrophe, or the 
short wave length instability, even though we have taken 
account of the surface tension effect by introducing the 
capillary length. 

Actually, numerical simulations of the equations even- 
tually yield a zigzag structure of the solution in the short- 
est discretization length; This cannot be accepted as a 
solution of the differential equations. 

Mathematically, the problem is that the surface ten- 
sion effect in ea. (|10|l is of the same form as k with the 
term that causes the instability in eq. lO , thus the surface 
tension effect never dominates to suppress the instability 
even in the short wave length limit, unlike in the case of 
MuUines-Sekerka/Saffman- Taylor problems. 

This may be seen in the linear stability analysis of the 
flat interface solution H13I) 11]. Since the flat interface 
solution is not steady, the linearized equations for the 
small deviation from it are not of constant coefficient, 
thus the analysis is not simple, but within the approxi- 
mation where the time variation is neglected during the 
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FIG. 2: Steadily propagating finger solutions without difltu- 
sion for R = 0.2 with V = 0.3, 0.2, and 0.1. (a) Shapes of 
the finger. The plots for diflerent speeds are shifted to avoid 
overlapping. Only the part near the tip is shown for each 
finger, (b) The line density of granules v.s. £. The finger tip 
is located aX £ = Q. 
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time scale of the perturbation (quasi-steady approxima- 
tion), the growth rate of the unstable perturbation 
with the wave number q in the unstable branch can be 
obtained as 



n{q) = 



1 



(21) 

where /j — f'{af{t)); Vf and af are given by ea. l(T^ . This 
is positive for any q > when /j < and 



n{q) 



R 



for q 



oo, 



(22) 



suggesting the ultra-violet catastrophe, thus the model is 
not well defined in the continuum limit. 



VI. DIFFUSION DRIVEN BY INTERFACE 
MOTION 

In order to avoid this instability, we introduce the dif- 
fusion term along the interface in Eq.© as 



da 
dt 



1 — KCT 



(23) 



with a new small length scale Hd- The adopted form 
represents the diffusion flux along the interface; The flux 
is proportional to the interface speed w„ . Such a diffusion 
flux is natural in the case where the diffusion is driven 
by the interface motion; The grains are driven randomly 
along the interface by the distance £d during the time 
(■D/vn, i-e. the time during which the interface moves by 
ti) in the normal direction. Note that we do not assume 
the diffusion perpendicular to the interface. 

The linear stability analysis within the same approxi- 
mation as above gives the perturbation growth rate 



m = \ 



-(^{R + £D)vtq^ - fi 



(113) for i? = 0.2 with = 0.2, 0.1, and 0.05. The steady 
solution of = 0.1 without the diffusion {£d — 0) is used 
as the initial state. Only the right halves of the inter- 
faces are shown and the time development is represented 
by the plots with the time interval At — 1. I have con- 
firmed the ultra-violet catastrophe is suppressed by the 
diffusion term; Except for the cusps at the edges of stick- 
ing regions, numerical solutions converges to a smooth 
solution as the smaller time and space discretization is 
used for integral in contrast with the diffusionless model, 
where the zigzag structure in the smallest discretization 
scale develops eventually all over the interface due to the 
ultra-violet catastrophe. In all cases of Fig|31 the fingers 
extend at a speed of 0.1 initially with keeping the initial 
shape as the steady solution does, but eventually the fin- 
ger becomes unstable and develops wavy structures. The 
length scale of the emerging structure is shorter for the 
smaller £ u , but they are much larger than £ o . We do not 
see any tendency that the steady solution becomes stable 
for small £d even though the initial configuration of the 
steady solution for the speed V = 0.1 with R = 0.2 is out- 
side the range Ea. lj^ . thus does not have a singularity 
in a at the tip. 

Fig^ shows the results for the case of i? = 0.2 and 
Id = 0.1. The initial states are the steady solutions of 
V = 0.1, 0.2, and 0.3 for £d = 0. One may see some dif- 
ferences in the way how the steady solutions are destabi- 
lized for different speeds, but again the steady solutions 
are unstable for all the cases, and complex developments 
of the boundary are seen after the instability. 

The unstable development seen in Figs. 3 and 4 shows 
some similarities to the pattern found in '3|; Concave 
parts of the interface evolve eventually into cusps with 
the tip size oi £d, and the interface are pinned by them. 
This leads to the irregular development of the interface 
although there is no randomness in the present model. 



(24) 



{R + £D)viq' -n) + 4 (-/fOtXf - £DRvfq^]viq^ 



which gives 



f2(g) < q> qs = 



therefore the catastrophe is suppressed. 




(25) 



VII. NUMERICAL SIMULATION OF FINGER 
SOLUTIONS 

I have performed simulations on the model with this 
diffusion in order to see if the finger solutions are stable. 

Let us start by examining the effect of the length scale 
£d- Figl^lshows the results of simulations of Eas. ljlOl) and 
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FIG. 3: The time developments of the interface for 7? = 0.2 
with £d = 0.2, 0.1, and 0.05. Only the right halves of the 
fingers are shown. The time sequences are shown with the 
time interval At = 1. The steady solution for V — 0.1 with 
£d = are used as the initial configurations for all the cases. 
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FIG. 4: The time developments of the interface for 7? = 0.2 
and Id = 0.1 from some initial configurations. The steady 
solutions of the speed V = 0.1, 0.2, and 0.3 with in = 
are used as the initial configurations. The time sequences are 
shown with the time interval At = 1. 

VIII. CONCLUDING REMARKS 

Before concluding, let us make some remarks. 

First of all, the sweeping phenomena can be seen com- 
monly; It is not limited to the specific experiment we re- 
ferred to, although there have not been many controlled 
experiments and theoretical analyses. Another example 
may be found in a pattern formation of deposit in a dry- 
ing droplet, 12]. Actual situations varies and may not be 
as simple as the one we have analyzed in this paper, but 
there should be a class of phenomena characterized by 
the sweeping phenomenon as is discussed here. 

Secondly, the boundary dynamics should be able to 
be derived from the phase field model as the narrow in- 
terface width limit. In the case of crystal growth and 
the viscous fingering, the interface width is the shortest 
length scale in the problem, and the boundary dynamics 
can be derived from the phase field model by taking ap- 



propriate limits [T3|. On the other hand, in the sweeping 
phenomenon, the length scale over which the granular 
density varies in the normal direction to the interface is 
in the same order with the interface width"?!, which com- 
plicates the formal derivation of the boundary dynamics. 

Lastly, some comments on the diffusion are in order. 
We have formulated the sweeping dynamics as the small 
diffusion limit, but in a real system, some form of diffu- 
sion should exist. We have also found that the diffusion- 
less model shows the short wave instability. In Ref.0], 
the ordinary diffusion is assumed within the interface, 
which results in patterns with a larger scale for a slower 
process; This does not seem to agree with the experiment. 
On the other hand, we regularized the model with the dif- 
fusion characterized by the short length scale £d, which 
may correspond to the grain size in the experiment^]. 
The diffusion along the interface is considered to be 
driven by the interface motion, thus the grains do not dif- 
fuse when the interface gets stuck, therefore, the slower 
interface motion does not lead to a larger scale. 

In summary, the boundary dynamics of the sweep- 
ing interface is constructed based upon geometrical and 
physical analysis of the process. To suppress the short 
length scale instability, the diffusion driven by the inter- 
face motion along the interface is introduced with the 
short length scale £d- In the case of no diffusion, we 
obtain the steadily propagating finger solutions for a fi- 
nite range of propagation speed, but the numerical sim- 
ulations with the small diffusion suggests that they are 
unstable, and a complex behavior of the interface is seen. 
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